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We discuss a solvable model describing an Anderson like impurity in a BCS superconductor. The model 
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Introduction 



The problem of magnetic impurities in superconductors is an old yet still very interesting problem. 
Starting with the work of Abrikosov and Gor'kov [1], numerous theoretical studies have addressed this 
problem. The models most frequently used to describe the magnetic impurities are the Kondo (s-d exchange) 
model [2] and some form of the Anderson model [3] . Most earlier works are based on a perturbative approach. 
A review of the earlier works and various references can be found in [4] . 

An important effect due to the impurities involves the existence of localized excited states within the 
superconducting energy gap. The results of Muller-Hartmann and Zittartz on gap states due to a Kondo 
impurity are well known [5] . More recent studies of the gap states include a numerical rcnormalization group 
analysis for a Kondo impurity [6], quantum Monte Carlo techniques applied to the symmetric Anderson 
model [7], large-N slave boson methods for an N-fold degenerate Anderson impurity [8] and a Bethe ansatz 
analysis for an Anderson impurity by Rupasov [9] . The closely related problem of an Anderson like impurity 
in a Hubbard chain, which has a spin gap, was studied in [10]. Recent experimental work on single impurities 
was reported in [11]. 

In this paper we consider a rather simplified version of the problem of a magnetic impurity in a BCS 
superconductor. The virtues of our model are that it is solvable by elementary methods, and also that it 
reveals some of the expected features of more complicated models. Our solution represents a self-consistent 
treatment of the BCS pair potential and a fully quantum impurity. We came to consider this model after a 
detailed study of the more difficult theory corresponding to a bulk mass term in the usual Kondo model [12]. 
Here, the interaction with the impurity is taken to be of the Anderson form. However unlike the Anderson 
model, the impurity 'd' orbital is considered to have no spin (or orbital) degeneracy, hence it is described by 
a single fermion operator d and its adjoint d\ with no spin index. Since the conduction electrons have spin, 
the impurity interaction for this model is intermediate to that of the Anderson model and the resonant level 
model [13] at the Toulouse point. We study the model by mapping it onto a boundary field theory (BFT) 
and solving for the reflection S-matrix (i.e. boundary scattering matrix). The BFT consists of two decoupled 
Ising models with boundary magnetic fields. The Ising fermions can be associated with the quasi-particles 
resulting from a Bogoliubov transformation. From this mapping one sees that whereas a Kondo impurity 
has no non-trivial flow in boundary conditions, an Anderson impurity does have a flow between free and 
fixed boundary conditions. A similar model, where the localized 'd' electrons do have a spin degeneracy, i.e. 
the Anderson model with no Coulomb repulsion (U — 0), was studied some time ago in [14] using the more 
conventional Green's functions approach. 

Our solution exhibits a single-impurity version of Andreev scattering. Andreev scattering occurs at 
the boundary between a superconductor and an ordinary metal, wherein quasi-particles are reflected into 
quasi-holes. For a single magnetic impurity, superconductivity can be destroyed locally in the vicinity of 
the impurity by the magnetic field of the impurity, in analogy with the Meissner effect. Thus one is led to 
expect that the impurity reflects quasi-particles into quasi-holes, and our solution confirms this. 

It is well-known that magnetic impurities compete with superconductivity in a manner analogous to 
the Meissner effect: if the concentration of impurities is large, enough bound states can open up in the gap 
that the gap disappears, destroying superconductivity. For a single impurity the competition is between the 
formation of Cooper pairs and the binding of electrons to the impurity, and is manifested in the properties 
of the bound states localized at the impurity. In our model, from the poles in the reflection amplitude for 
each Bogoliubov quasi-particle we can determine the bound states of the quasi-particle with the impurity. 
Simple energetics determines whether the bound state is stable. The impurity has a magnetic energy scale 
Eimp- When the gap A is zero, this energy scale appears as a resonance in the reflection S-matrix, and 
corresponds to an unstable bound state of the original electrons with the impurity. (In the doubled boundary 
Ising description, this resonance corresponds to a splitting in the degeneracy of the two Ising ground states 
induced by the boundary magnetic field.) When the gap is non-zero, the state of energy Ei mp can become 
stable if it is forbidden to decay into the bulk superconductor because of the gap; this occurs if A > E lmv . 
This bound state has an excitation energy, Eb, less than the superconducting gap A. The bound state energy 
is a function of E- lnip and A, and is given in (3.15) below. We find that the bound states in the gap disappear 
when the strength of the impurity is large enough, i.e. when E lmv > A. In this case the 'Meissner effect' is 
reversed and superconductivity is favored: the quasi-particle-impurity bound state becomes unstable and it 
is more energetically favorable to form Cooper pairs. 
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1. The Model 



The model we use to describe an impurity in a superconductor consists of the (mean-field) BCS Hamil- 
tonian and an Anderson like interaction, 

H = H BCS +H I (1.1) 
^bcs = £ ^ C£r - A £ (4 t cl Sl + c_^) (1.2) 

fecr 

where ct are conduction electron creation operators of momentum k and spin a (f up and J. down), and A 

is the energy gap. The impurity is at the origin in three spacial dimensions. The operators d and d) describe 
a spinless fermion localized on the impurity, satisfying 

{d,cft} = l, rf 2 = dt 2 = . (i.4) 

These operators act on a two-dimensional Hilbert space |0) d , |l) d where |l) d = d j( \0) d . In the unperturbed 
theory these two states are taken to be degenerate in energy. Since the impurity level has no spin degeneracy, 
there is no Coulomb repulsion term as in the Anderson model. The energy £g is measured relative to the 
fermi energy e F , £g = eg — e F , where eg is the electron kinetic energy. In the impurity interaction Hj, Vg is 
the hybridization matrix element between the band states and the impurity states. 
We can diagonalize the BCS part by a Bogoliubov unitary transformation [15] 

Ctt = Urbr. + aT-vrbt , C 1 -. = urbt — a-rvrbt, , (1.5) 
fcj fc fel fe k k2 — fej. k k2 k fe fel' v 1 

with 



'fe 



l-*<lfl-M-fr, Jjg>£ (i.T) 

and 

^=^1 + A2 - (L8) 
In terms of the quasi-particle operators (frt . , br . ) , which are often referred to as Bogoliubons, -ffecs becomes 

KJ J 

ffBCS=E¥k g fc' (L9) 

fej 



where 

C -wg, if |fc| < k F - (1 1Q) 



[WjJ, if|fc|>fc F . 

The ground state of (1.9), |0) 6 , is obtained by filling all levels with \k\ < k F . Rather than working with the 
b s, it is more convenient to introduce the following 'modified' quasi-particle operators 



b = \hv ti\k\>k F \ b f&fe, if\k\>k F ; 
fel Hg 2 , if|fc|<fc F , fe2 \b Uv if|fc|<fc F . 



(1.11) 



These operators satisfy canonical commutation relations and will simply be referred to as the quasi-particles 
of the system. The BCS Hamiltonian maintains its diagonal form (1.9) when expressed in this new basis. 



The Bogoliubov transformation applied to the interaction gives 

Hl = ^ E Hi b ki)( d + dt ) + ( fo L + b k2)( d rft )] . (1-12) 

k 

where we have set k — kp in Vg, = Vk F , and in the coefficients of the unitary transformation (1.6). In 
(1.12) we have kept the terms bt and &g d, which were omitted in the 'rotating wave approximation' made 
in [9]. 

Assuming the impurity to be point-like with only s-wave scattering, we expand the operators (£>t , fer.) 

KJ J 

in spherical harmonics and retain only the angular momentum I = m = terms (the / ^ terms do not 
couple to the impurity). This gives the one dimensional Hamiltonian (sum over j) 

H = F CC ^ {u{k)b){k)b 3 {k) + X= [(b\(k) h(k))(d + dt) + (fet(fc) + b 2 {k)){d dt)] ) , (1.13) 



where 



We have linearized £g about the fermi surface, £g = |/c| — &f, setting vf = 1, and defined k = \k\ — kp to be 
the momentum relative to the fermi surface. The operators &j(fc) create quasi-particles of momentum k + kp 
and satisfy 

{bi(k),b](k')} = 2ir5ij5(k — k'). (1.15) 

Note that V has dimension ^/energy; below we show that the scale E lmp introduced in the introduction is 
simply V 2 . 

A final transformation ib\(k) — > b\(k) and — ib\(k) — > &j(fc), allows iJ to be rewritten as 

H = H 1 +H 2 (1.16) 

^ = (w(A)65(fc)6 j (A ! )+i^(&5(fc)+6 J -(A ; ))a J -) (j = l,2), (1.17) 

where 

d = (d + <**), 02 = -*(<*-<**). (1.18) 
The operators dj anticommute with the 6 operators and satisfy 

a] = aj, a| = 1, {ai,a2} = 0, (1-19) 
implying if 2 ] = 0. Thus we only need to concentrate on one copy and henceforth drop the subscript. 

2. Mapping onto a Boundary Field Theory 

In this section we map (1.17) onto a boundary Ising field theory. As in the usual Kondo model, s-wave 
scattering has reduced the problem to a one-dimensional problem on the half-line r > 0, where r is the 
spherical radial coordinate. Henceforth, r is designated as the spacial variable —x. 

The action for an Ising field theory on the half-line with a magnetic field at the boundary is [16] 

S = Sbulk + S'a-frcc + Sbc + Sint (2.1) 

i r°° r° _ _ _ 

Sbuik = 2 / dt dx 0(#t - W + yj(d t + d x )4> + 2A^) (2.2) 

^ J — oo J — oo 
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if 00 — i f 00 

S hc = -- dtiH>\ x=0 , S a - bcc = - dta(t)d t a(t) (2.3) 

V f°° — 

S mt = -i-j=J dt(^ + iP)a(t)\ x=0 . (2.4) 

The term Sbuik is simply a free Majorana action with A being the fermion mass, »!?bc serves to enforce the 
free boundary condition 

■ip^Tp atar = 0, (2.5) 

for V = 0, and Si n t gives the interaction between the boundary spins and the boundary magnetic field V. 
The fermionic field a(t), which anticommutes with (ip 7 ip), describes the ground state degeneracy due to the 
two different expectation values of the spin field and satisfies a 2 = 1. Thus 5 a -frcc is the kinetic term for 
a(t). The bulk equations of motion are 

(d t - d x )^ + a^ = o, (d t + d x y$ - = o. (2.6) 

Mode expansions satisfying (2.6) can be written as 



(2.7) 



where 

{A{k),AHk')} = 2nS(k-k') (2.8) 



and k ■ x — uj{k)t — kx, cj(k) = \Jk 2 + A 2 . 

The Ising Hamiltonian associated with (2.1) takes the form 

H = r S ( u( ^ tWi(fc) + ^ ( g {k)A{k) + 9*(k)A\k)) a)j , (2.10) 

where 

g(k) = u(k)+v(k). (2.11) 

The second term in (2.10) gives the boundary interaction. Comparing (2.10) with (1.17), we can map each 
Hj onto the Ising Hamiltonian if we relate the b operators to the A operators as follows 

b(k) = g{k)A{k) if k > (2.12a) 

b(k) = g*{k)A\k) if k < 0, (2.12b) 

and identify a,j with a. The quasi-particles correspond to the Ising fermions and the impurity coupling V 
plays the role of a boundary magnetic field. Using 

g(k)g*(k) = 1, (2.13) 

it is easily shown that the b operators satisfy canonical commutation relations 

{b{k),b\k')} = 2n5{k-k'). (2.14) 

3. Quasi-Particle Bound States 

Having mapped (1.17) onto the BFT (2.1), the impurity interaction can then be described by the 
reflection S-matrix (boundary scattering matrix) for the fermion operators. We compute the reflection 
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amplitude from the action as discussed in [16]. The boundary terms in the variation of (2.1) lead to the 
boundary equation of motion (at x = 0) 

-8 t {^-^) = V 2 (ip+^). (3.1) 

Here we see the flow between free (V = 0) and fixed (V = oo) boundary conditions. Interpreting (3.1) 
to hold when acting on a formal boundary operator B, we substitute the mode expansions and define the 
reflection S-matrix, or reflection amplitude, R(k), through 

A\k)B = R{k)A\-k)B. (3.2) 

The reflection amplitude for the Ising model was computed in [16]. Let us parametrize the energy and 
momentum with the rapidity 9 

uj{9) = Acosh(9, k{6) = Asinh<9. (3.3) 



Define A{9) = A{k)^/U{k), satisfying {A(6), A^{9')} = 2nS{9 - 6'). We find 

A\6)B = R(9)AH-9)B, (3.4) 

where 

The amplitude R b is the reflection amplitude for the Bogoliubov operators b(k) 

tf(k>0)B = R b (9)b(-k)B, b(k < 0)B = R b (9p (-k)B. (3.6) 

Defining the quasi-electron operators e(k) for k > and quasi-hole operators h(k) for k < 

e(k) = b(k), forfc>0; h(k) = tf(k), for k < 0, (3.7) 

then (3.6) can be expressed as 

e\k)B = R b {k)h)(-k)B, h\k)B = R b {k)e\-k)B . (3.8) 

The quasi-electrons scatter into quasi-holes and vice-versa at the boundary. This is a form of Andreev 
reflection, which usually occurs at the boundary between a superconductor and a normal metal. Here the 
boundary is a single impenetrable impurity, i.e. there is no transmission. For the original quasi-particles, 
the 6's, equation (3.8) implies that quasi-electrons of one type, say b\(k), scatter into quasi-holes of the other 
type b 2 {-k). 

In terms of the energy E(9) = u)(6), R b can be written as 

E 2 + (V 2 ) 2 - 2AV 2 

Rb(E) = { ' — (3.9) 

(V 2 ) 2 -E 2 - 2(V 2 - A)(A + i^E 2 - A 2 ) 

For A = 0, (3.9) becomes 

K(E) = - E^-i-iTE ' <""> 



where 



E hap = V, T = 2V Z . (3.11) 



The energy E- lmp is a 'magnetic energy' scale associated with the magnetic impurity. This becomes clearer in 
the Ising picture where V effectively acts as a local magnetic field. Near E w E- lmp , R b {E) takes a Lorentzian 
form 
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From (3.12) we see that when A = 0, there is a resonance in the reflection amplitude at E = Ei mp with a 
width r. This resonance corresponds to an unstable electron- impurity state with an energy Ei mp and decay 
rate T. If the gap A is non-zero, then this state can become a stable bound state (see below). 

As with the bulk S-matrix, poles of the reflection S-matrix provide information on bound states. Poles 
in the physical strip < Im# < 7r can be indicative of either bulk bound states or boundary bound states. 
Since the bulk theory is a free fermion theory with trivial scattering (S = —1), we only have to consider the 
latter case. The quasi-particle reflection amplitude is 

_ »sinhg + (£ imp /A-l) 
Rb{9) zsinh0-(iWA-l)- ( } 

We will now express all quantities in terms of the physical energy scale E lmv . Writing as iu, poles of Rb(0) 
occur when 

sinu = l_^£. (3.14) 

Restricting u to be in the physical strip, there is a single pole at u — u b satisfying (3.14) provided Ei mp < A. 
(Note that there is also a pole at tib = tt — Ub, corresponding to scattering in the cross channel. The results 
that follow also apply to this pole.) Thus a quasi-particle forms a bound state with the impurity if Ei mp < A. 
This is the same state that appears as a resonance in (3.12) for a vanishing gap; it has become stable because 
its decay into the bulk is forbidden if E; mp < A. The excitation energy (above the ground state) of the 
bound state, Eb, is 

E b = A cos u b = ^ AE imp (2 - . (3.15) 

For £i m p < A, the excitation energy lies within the superconducting gap (Eb < A) and the bound state is 
stable. When Ei mp > A, the excitation energy falls in the continuum and there is no bound state. 

If i?i m p = 0, the pole occurs at itt/2, implying that the ground state is degenerate. This makes sense 
since for zero interaction the ground states are 

|0> 6 ®|0> d and |0) 6 ®|l) d , (3.16) 

where |0) b is the quasi-particle ground state. In terms of the quasi-electron and quasi-hole operators 
(e(k),h(—k)) with k > 0, \0} b is a zero energy state satisfying 

e(fc)|0> 6 = M-*)|0> 6 = 0. (3-17) 

As Eimp increases above zero, the degeneracy is lifted by the interaction. The states (3.16) are split into 
a ground state and the bound state with energy (3.15). In the Ising picture the degenerate ground states, 
|0, ±) 7 , are labeled by the expectation values of the spin field a(x), (<r(x)}± = zta, where a is the spontaneous 
magnetization. A non-vanishing magnetic field, ^/E; mp > 0, removes the degeneracy of the free boundary 
spin states |0, ±) 7 . 

Based on the above discussion, we can interpret Andreev reflection of the quasi-particles as follows. 
Suppose we extend our system by adjoining a lattice of impurities for x > 0. In the region x > supercon- 
ductivity is more or less destroyed and we have a superconductor-normal conductor boundary at x = 0. In 
this case there will be Andreev reflection at the boundary, as well as transmission. A lattice consisting of a 
single impurity then behaves as a superconductor-normal conductor barrier with the probability for Andreev 
reflection unity, |i?;,| 2 = 1, just as would be the case for subgap particles incident on a superconductor. 

So far we have only been dealing with one of the Hamiltonians Hj in (1.16). Each type of quasi-particle, 
b\(k) or b\(k), will form a bound state. If we represent the ground state of Hj by |0) ■ (for small Ei mp > 0) 
and the bound state by |£>) ., then the ground state of H = H\ + H 2 , |0), will be 

|0) = |0) 1 ®|0) 2 , (3.18) 

and the excited states are 

|1) = |6) 1 ®|0) 2 , |1'> = |0> 1 ®|6) 2 , (3.19) 
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both with the same energy and lying inside the gap. 



4. Conclusions 

We have shown that the model (1.1) is a solvable system, which is equivalent to two decoupled Ising 
models with boundary magnetic fields. From the boundary field theory we calculated the reflection am- 
plitudes for the quasi-particles. At the boundary the quasi-particles are Andreev reflected. Provided the 
magnetic energy E- imp (or equivalently the hybridization matrix element/magnetic field V) is not too large, 
bound states exist, consisting of a quasi-particle bound to the impurity. The energy of the stable bound 
states is less than the gap, Eb < A. 

We hope this calculation further motivates the use of techniques common to integrable systems for the 
study of impurities in superconductors. It would be interesting to apply the techniques discussed here, and 
in [12], to more complicated systems, such as the original Anderson model, or magnetic impurities in d-wave 
superconductors studied in [17]. 
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